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ABSTRACT 



Aims. Our aim is to search for and study the origin of the low-amplitude and long-periodic radial velocity (RV) variations in K giants. 
Methods. We present high-resolution RV measurements of K2 giant e CrB from February 2005 to lanuary 2012 using the fiber-fed 
Bohyunsan Observatory Echelle Spectrograph (BOES) at the Bohyunsan Optical Astronomy Observatory (BOAO). 
Results. We find that the RV measurements for e CrB exhibit a periodic variation of 417.9 + 0.5 days with a semi-amplitude of 129.4 
± 2.0 m s . There is no correlation between RV measurements and chromospheric activity in the Ca II H region, the Hipparcos 
photometry, or bisector velocity span. 

Conclusions. Keplerian motion is the most likely explanation, with the RV variations arising from an orbital motion. Assuming a 
possible stellar mass of 1.7 ±0.1 M Q for e CrB, we obtain a minimum mass for the planetary companion of 6.7 + 0.3 M Sup with an 
orbital semi-major axis of 1.3 AU and eccentricity of 0.11. We also discuss the implications of our observations for stellar metallicity 
versus planet occurrence rate and stellar mass versus planetary mass relations. 
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1. Introduction 

As a star evolves towards the red giant branch, its rotational ve- 
locity decreases. Simultaneously, the stellar surface begins to 
cool significantly, and it may develop surface inhomogeneities 
such as stellar spots. The activity can cause changes in spectral 
line profiles, which develop in a complex manner with the grav- 
itational influence of a planetary companion. Thus, this and stel- 
lar pulsations make it more difficult to detect exoplanets around 
giants than around main-sequence (MS) stars. Nevertheless, cool 
evolved stars are suitable samples for precise RV measurements 
and for studying exoplanets because they offer new clues to help 
develop theories of planet formation and help explain the long- 
term survival of planetary systems. 

Since the discovery of the first planetary companion around 
K giant star l Dra (Frink et al. 2002), about 30 planetary compan- 
ions have been detected so far around evolved K giant stars using 
the precise RV method (Setiawan et al. 2003a, 2003b; Mitchell et 
al. 2003; Hatzes et al. 2005, 2006; Reffert et al. 2006; Ddllinger 
et al. 2007, 2009; Johnson et al. 2008, 2010, 2011; Sato et al. 
2008a, 2008b, 2010; de Medeiros et al. 2009; Niedzielski et al. 
2009; Han et al. 2010; Lee et al. 2011; Moutou et al. 2011; 
Wittenmyer et al. 2011; and Gettel et al. 2012). They seem to 
show different properties than MS stars. Giants that host plan- 
ets do not favor high-metallicity objects (Pasquini et al. 2007) 
and trend to harbor more massive planets (Lovis & Mayor 2007; 
Dollinger et al. 2009). 

For the past seven years, we have conducted precise RV mea- 
surements of 55 K giants. Here, we present a long-period and 
low-amplitude RV variation of K giant e CrB. In Sect. 2, we 
describe the properties of e CrB and analyze the periods of RV 



measurements in Sect. 3. Possible origins of the RV variations 
are investigated in Sect. 4. Last, we discuss the implications of 
this work in Sect. 5. 



2. The properties of e CrB 

The cool evolved star of spectral type of K2 III (ESA 1997), e 
CrB (= HD 143107 = HR 5947 = HIP 78159) has an apparent 
magnitude of 4.14. The main parameters come from Massarotti 
et al. (2008), who present rotational velocities for a sample of 
761 giants selected from the Hipparcos catalog. They measured 
new rotational velocities for all the sample using spectroscopic 
line broadening methods with the CfA Digital Speedometers 
(Latham 1992). They estimated the rotational velocity of e CrB 
to be 2.4 km s . Based on the rotational velocity of 2.4 km s _1 
and stellar radius of 21 R (Massarotti et al. 2008), we derived 
the range for the estimated velocity of the rotational period as 
P mt = 2nRi, I (v rot sin i) = 442.7 days. 

The atmospheric parameters of e CrB were determined based on 
the program TGVIT (Takeda et al. 2005). We used 245 equiv- 
alent width (EW) measurements of Fe I and Fe II lines, which 
resulted in T eff = 4406 ± 15 K, log g = 1.94 ± 0.08, v micl0 = 1.59 
+ 0.07 km s _1 , and [Fe/H] = -0.094 ± 0.001. 

The stellar mass has been determined somewhat differently 
by individual authors (because it is very sensitive to log g). We 
estimated the stellar mass from the theoretical stellar isochrones 
using its position in the color-magnitude diagram based on 
Bertelli et al. (1994) and Girardi et al. (2000). We adopted 
a version of the Bayesian estimation method by J0rgensen & 
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Table 1. Stellar parameters for e CrB. 



Parameter 


Value 


Reference 


Spectral type 


K2III 


Hipparcos 


m v [mag] 


4.14 


Hipparcos 


B-K[mag] 


1.23 


Hipparcos 


age [Gyr] 


1.74 ± 0.37° 


Derived 


Distance [pc] 


67.9 


Shaya& Oiling (2011) 


RV [km s" 1 ] 


- 30.92 


Massarotti et al. (2008) 


Parallax [mas] 


13.97 ± 0.68 


Kharchenko et al. (2009) 


Diameter [mas] 


3.08 ± 0.16* 


Richichi & Percheron (2002) 


Teff [K] 


4365 ± 9 


Massarotti et al. (2008) 




4406 ± 15 


This work 


[Fe/H] 


-0.32 


Massarotti et al. (2008) 




-0.094 ± 0.001 


This work 


logg 


2.3 


Massarotti et al. (2008) 




1.94 ± 0.08 


This work 


R* [Re] 


21 


Massarotti et al. (2008) 


M* We] 


1.7 + 0.1" 


Derived 


U [L Q ] 


151 


Massarotti et al. (2008) 


v 10t sin i [km s J ] 


2.4 


Massarotti et al. (2008) 


P rot / sin ( [days] 


442.7 


Derived 


Vraicro [km S~'] 


1.59 ± 0.07 


This work 



<a> Derived using the online tool (http://stevoapd.inaf.it/cgi-bin/paraml. 
(i) Computed on the basis of a formula as described in method, and 
error is included whenever possible. The uniform disk angular diameter 
is 2.73 ± 0.03 mas, and the limb-darkened angular diameter is 2.80 ± 
0.03 mas. 



Lindegren (2005) and da Silva et al. (2006) using the determined 
values for r e ff, [Fe/H], M v , and parallax. Our estimated param- 
eters yield a stellar mass of 1.7 + 0.1 M . The basic stellar pa- 
rameters of e CrB are summarized in Table [TJ 

3. Observations and analysis 

Since 2003, we have conducted a precise RV survey of 55 
K giants using the fiber-fed high-resolution (R = 90 000) the 
Bohyunsan Observatory Echelle Spectrograph (BOES; Kim et 
al. 2007) attached to the 1.8-m telescope at Bohyunsan Optical 
Astronomy Observatory (BOAO) in Korea. 

We acquired 52 spectra for e CrB from February 2005 to 
January 2012. Each estimated signal-to-noise ratio (S/N) at an 
iodine (I2) wavelength region (4900 - 6000 A) is about 250 with 
a typical exposure time ranging from 240 to 480 seconds. An I2 
absorption cell was used to provide precise RV measurements. 
The extraction of normalized 1-D spectra was carried out with 
IRAF (Tody 1986) software and DECH (Galazutdinov 1992) 
code. The I2 analyses and precise RV measurements were un- 
dertaken using the RVI2CELL (Han et al. 2007) developed at 
the Korea Astronomy & Space Science Institute (KASI). The 
RV measurements for e CrB are listed in Table [2] 

The RV standard star r Ceti shows a stable RV with an rms 
scatter of 6.8 m s _1 (Figure [TJ over the time span of our obser- 
vations and demonstrates the long-term stability of the BOES 
(Lee et al. 2011). RV measurements of e CrB show a standard 
deviation of 76.4 m s , which is 1 1 times more than for the RV 
standard star (Figure [TJ. 

Period analysis is applied in order to determine the signifi- 
cance of periodic trends or variabilities in the time series. The 
Lomb-Scargle periodogram (Lomb 1976; Scargle 1982), a use- 
ful tool for investigating long-period variations for unequally 
spaced data, was calculated for the RV time series for e CrB. 



Table 2. RV measurements for e CrB from February 2005 to 
January 2012 using the BOES. 



JD 


ARV 


+<T 


JD 


ARV 


+<T 


-2 450 000 


ms~' 


m s 


-2 450 000 


m s 


m s 


3430.324001 


-81.2 


6.7 


4994.182859 


-48.5 


17.8 


3460.235768 


2.5 


7.4 


4995.208262 


-45.3 


6.9 


3460.297588 


-9.7 


7.0 


4996.099283 


-36.8 


7.0 


3506.253786 


130.3 


8.2 


4996.106632 


-23.7 


9.4 


3545.042399 


123.5 


7.3 


5130.932411 


-4.1 


8.3 


3730.408100 


-65.3 


8.1 


5248.368418 


173.7 


6.6 


3818.237686 


-68.4 


6.2 


5321.137848 


144.6 


6.8 


3821.138258 


-80.1 


7.3 


5356.148060 


48.7 


7.9 


3891.145666 


43.3 


8.6 


5810.927302 


50.3 


7.1 


3899.184958 


60.4 


7.3 


5810.993537 


27.4 


6.7 


4147.298612 


-57.9 


6.7 


5820.022432 


-61.3 


6.8 


4210.246254 


-83.0 


9.2 


5820.028150 


-62.0 


7.3 


4262.234769 


-100.3 


8.7 


5820.034017 


-64.2 


8.0 


4264.223335 


-25.9 


7.3 


5841.950835 


-53.2 


6.6 


4458.409299 


150.6 


6.9 


5841.955858 


-52.4 


6.7 


4483.406189 


158.4 


8.7 


5841.960546 


-50.3 


6.6 


4505.293538 


113.9 


8.0 


5842.942221 


-54.5 


7.0 


4516.321792 


86.8 


7.8 


5842.942221 


-54.5 


7.0 


4536.294962 


17.9 


7.5 


5842.942221 


-54.5 


7.0 


4538.204751 


33.1 


7.8 


5911.396383 


-33.0 


16.8 


4619.213839 


-138.5 


18.8 


5911.401059 


-25.3 


20.0 


4847.360348 


149.9 


8.0 


5911.406117 


-50.1 


18.8 


4930.186145 


83.3 


6.9 


5933.351526 


-23.0 


8.8 


4970.217174 


-5.5 


8.5 


5933.355508 


-17.1 


8.9 


4970.224605 


-10.7 


8.1 


5933.359501 


-26.2 


7.9 


4971.157049 


-6.9 


7.3 


5933.363483 


-26.9 


7.7 



Figure [2] (top panel) shows a significant power at f\ = 0.00239 
cd" 1 (P = 417.9 days). 

A false alarm probability (FAP) is a metric that expresses 
the significance of a period. A false alarm arises in period anal- 
ysis techniques when a period is incorrectly found where none 
exists in reality. We determined the significance of a period by 
calculating an FAP for the dominant period by a bootstrap ran- 
domization technique (Press et al. 1992; Kiirster et al. 1999). 
We found an FAP of less than 10~ 5 for the period of 418 days. 
Figure [2] (bottom panel) shows the Lomb-Scargle periodogram 
of the residual after subtracting the frequency of f\ = 0.00239 
c d _I , and we were able to find a remaining peak at the frequency 
of /1 = 0.00497 c d" 1 (P ~ 201 days) with an FAP of 1.4 + 0.2%. 

4. Origin of the RV variations 

RV of giant stars may be affected by surface phenomena and by 
the presence of an orbiting planet. While short-term RV vari- 
ations (hours to days) with a low-amplitude in evolved stars 
have been known to be the result of stellar pulsations (Hatzes 
& Cochran 1998; Hekker et al. 2006), long-term RV variations 
(hundreds of days) may be caused by three kinds of phenomena: 
stellar pulsations, rotational modulations by inhomogeneous sur- 
face features, or planetary companions. 

4. 1 . Orbital mode 

The source e CrB reveals a periodic signal in our RV variations, 
and we find the best-fit Keplerian orbit with a P = 417.9 + 0.5 
days, a semi-amplitude K = 129.4 + 2.0 m s , and an eccentric- 
ity e — 0.11 + 0.03. Figure [TJ shows the RV curve as a function 
of time for e CrB and the residual after extracting the main fre- 
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Fig. 1. Variations in the RV standard t Ceti (to/? panel), RV curve 
(middle panel), and rms scatter of the residual (bottom panel) for 
e CrB from February 2005 to January 2012. The solid line is the 
orbital solution with a period of 418 days and an eccentricity of 
0.11. 
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Fig. 2. The Lomb-Scargle periodogram of the RV measurements 
for e CrB. The periodogram shows a significant power at a fre- 
quency of 0.00239 c d -1 corresponding to a period of 417.9 days 
(top panel). And after subtracting the main frequency variations 
(bottom panel), the highest significant power at a frequency is 
0.00497 c d" 1 (P = 201.2 days). The horizontal dashed line in- 
dicates an FAP threshold of 1%. 



quency. As can be seen in Figure [TJ the rms of the RV residuals 
are 25 m s -1 , which is significantly larger than the RV precision 
of 6.8 m s and than the typical internal error in individual e 
CrB RV accuracy (~ 8 m s _1 ). There are no obvious peaks in the 
residual periodogram (Figure|2]i. It is a general tendency that the 
rms of RV of a median value of 20 m s _1 is shown in K giants 
(Hekker et al. 2006). We derive the minimum mass of a plane- 
tary companion m sin i = 6.7 ± 0.3 Mj up at a distance of a = 1.3 
AU from e CrB. All the orbital elements are listed in Table[3] 



4.2. Chromospheric activity 

The EW variations of Ca II H & K lines are sensitive to atmo- 
spheric activity and are frequently used as chromospheric activ- 
ity indicators. The emissions in the Ca II H & K core are formed 
in the chromosphere and show a typical central reversal in the 
chromospheric activity (Pasquini et al. 1988; Saar & Donahue 
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Fig. 3. Ca II H spectral region for e CrB. The dotted lines are 
individual spectra, and the red solid line is the mean spectrum. 



1997). The existence of extra emission at the center of the line 
implies that the source function in the chromosphere is greater 
than in the photosphere. Unfortunately, because the e CrB data 
do not have a sufficient S/N in the Ca II K line region, we used 
the Ca II H line region, which has a bit higher S/N than does the 
Ca II K line region. Although the S/N is insufficient, the Ca II H 
line is enough to estimate the emission feature in the line core, 
and it does not show any visible emission features. Among the 
52 samples of Ca II H line, 35 spectra with relatively high S/N 
were drown for comparison. Figure [3] shows individual features 
in the Ca II H line region, and the red line is the mean spectrum. 
Although there is no prominent emission in the central region, it 
exhibits some scatter, which could be from either noise or chro- 
mospheric variations. 

We also measured the variations in the Ca II H line EW, 
which was estimated in a spectral range of 3967.7 - 3969.3 A 
centered on the core of the Ca II H line to avoid nearby blend- 
ing lines. The mean EW is measured to be 1464.0 + 27.8 mA. 
The rms of 27.8 mA corresponds to a 1.9% variation in the EW. 
Figure |4] shows Ca II H EW variations as a function of time and 
the Lomb-Scargle periodogram. There is significant power at the 
frequency of 0.00028 c d -1 , corresponding to a period of 3579 
days, which is unrelated to the RV period of 418 days. 

4.3. Hipparcos photometry 

We analyzed the photometry to search for possible brightness 
variations due to the rotational modulation of cool stellar spots 
or due to the stellar pulsations. The photometry database avail- 
able for e CrB comes from 180 Hipparcos measurements from 
November 1989 to March 1993. The data maintained a photo- 
metric stability down to an rms scatter of 0.0066 magnitude cor- 
responding to 0.15% variations over the time span of the ob- 
servations. Figure|5]shows the Hipparcos photometric variations 
and the Lomb-Scargle periodogram of the measurements. There 
are no significant peaks near the period of 418 days. 

4.4. Line bisector variations 

Stellar rotational modulations by surface features can create 
variable asymmetries in the spectral line profiles, as well as RV 
variations. Thus, the variations in the shapes of spectral lines 
may identify the origin of RV variations. The difference in the 
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Fig. 4. The examinations of the origin of the chromospheric ac- 
tivity (Ca II H line) for e CrB. It shows JD vs. Ca II H EW 
variations (top panel) and the Lomb-Scargle priodogram of the 
variations (bottom panel). The vertical dashed line marks the lo- 
cation of the periods of 418 days. 



m 



6000 

JD-2450000 




0.00 0.01 0.02 0.03 

Frequency (c/d) 



Fig. 6. The examinations of the origin of the BVS variations for 
e CrB. It shows JD vs. BVS variations (top panel) from February 
2005 to January 2012 and the Lomb-Scargle priodogram of the 
BVS variations (bottom panel). The vertical dashed line marks 
the location of the periods of 418 days. 




Table 3. Orbital parameters of the RV measurements for e CrB. 



8200 8400 8600 

JD-2440000 
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Frequency (c/d) 

Fig. 5. The examinations of the origin of the Hipparcos photo- 
metric variations for e CrB. It shows the Hipparcos photomet- 
ric measurements (top panel) from November 1989 to March 
1993 and the Lomb-Scargle priodogram of the variations (bot- 
tom panel). The vertical dashed line marks the location of the 
periods of 418 days. 



width of the line at the top and at the bottom of the line profile is 
defined as bisector velocity span (BVS). 

We measured the BVS using the least-squares deconvolution 
(LSD) technique (Donati et al. 1997; Reiners & Royer 2004; 
Glazunova et al. 2008), which is calculated by the mean pro- 
file of the spectral lines. We also used the Vienna Atomic Line 
Database (VALD; Piskunov et al. 1995) to prepare the list of 
comparison lines. A total of ~ 3000 lines within the wavelength 
region of 4500 - 4900 A were used to construct the LSD pro- 
file. It excluded spectral regions around the I2 absorption region 
(4900 - 6000 A), hydrogen lines, and regions with strong con- 
tamination by terrestrial atmospheric lines. Finally, we estimated 
the BVS of the mean profile between two different flux levels of 
central depth levels of 0.8 and 0.25 as the span points. BVS vari- 
ations as a function of time and the Lomb-Scargle periodogram 
of the BVS for e CrB are shown in Figure[6] There are no obvi- 
ous periodic variations in the BVS measurements. 



Parameter 


Value 


Period [days] 


417.9 ± 0.5 


T^permstron [JD] 


2451235.3 ± 39.7 (1999.152) 


K [m s" 1 ] 


129.4 ± 2.0 


e 


0.11 ±0.03 


co [deg] 


133.1 ±33.3 


f(m) [M ] 


(9.227) x 10~ 8 


a sin i [AU] 


(4.943) x 10- 3 


cr(O-C) [m s" 1 ] 


25.0 


withM* = 1.7 ±0.1 [M ] 




m sin i [Mjup] 


6.7 ± 0.3 


a [AU] 


1.3 



5. Discussion 

From the analysis of the seven-year precise RV measurements 
for e CrB, we found strongly periodic variability in RV mea- 
surements with a period of 418 days and a semi-amplitude of 
129 m s _1 . There is no correlation between the RV measure- 
ments and the chromospheric activity, the Hipparcos photome- 
try, or the BVS. Thus, RV variations caused by a planetary com- 
panion seem the most likely explanation for the observed RV 
variations. 

We considered two interesting questions related to planets 
in evolved stars. First, a stellar chemical composition appears 
to be a major indicator of its probability of hosting planets. MS 
stars hosting planets are metal-rich (Santos et al. 2004; Fischer 
& Valenti 2005), and a survey of 160 metal -poor MS stars give 
no evidence of exoplanets (Sozzetti et al. 2009). In contrast to 
MS stars, giants hosting planets follow the same age-metallicity 
distribution as giants without planets (Pasquini et al. 2007). 
Recently, unlike the claim of Pasquini et al. (2007), Quirrenbach 
et al. (201 1) suggest that there is a clear trend for the planetary 
companion frequency to increase with metallicity just as with 
MS stars. However, since the sample selection criteria might 
have a significant influence on the correlations between com- 
panion frequency and metallicity in addition to the not having 
enough giants known so far, it would be too premature to reach 
a conclusion. 
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Fig. 7. Average mass of planetary companions as a function of 
stellar mass, taking all planets known as of February 2012 into 
account. 

Second, although the sample of intermediate-mass stars 
hosting planets is still limited, there could be a general tendency 
between the stellar mass and planetary companions being dis- 
played. It gives useful information on the yield of the planetary 
formation process. We consider six bins of stellar masses with 
equally spaced stellar mass of 0.5, up to 3 M , above which 
samples are too scarce, and compute an average mass of plane- 
tary companions in each bin. Figure [7] shows the results in a his- 
togram, where more massive stars do form more massive plan- 
etary companions (Johnson et al. 2007; Lovis & Mayor 2007; 
Dollinger et al. 2009), and an average planetary companion mass 
clearly increases in the range of stars with more than 2.5 stellar 
mass. Massive planetary or brown-dwarf mass objects originate 
in the circumstellar disk, as for giant exoplanets (Niedzielski et 
al. 2009). In principle, such a system could form from a suf- 
ficiently massive protoplanetary disk by means of the standard 
core-accretion mechanism (Ida & Lin 2005). Because more mas- 
sive stars probably have more massive disks and higher surface 
densities, larger amounts of matter can be accreted; thus, the fre- 
quency of massive planets could be higher. With 1 .7 M and a 
companion of 6.7 Mj up , e CrB seems to be consistent with this 
speculation. However, there is the obvious observational bias 
that it is more difficult to detect low-mass planetary compan- 
ions around a more massive star. Therefore, confirmation of this 
suggestion requires more discoveries of systems like e CrB. 

There has been an attempt to link larger disk masses of more 
massive stars to a possibility that the frequency of planets around 
giants is less dependent on metallicity than is the case for dwarf 
stars (Ghezzi et al. 2010). In that work, the lower metallicity for 
the giants hosting planets in comparison to the dwarfs hosting 
planets is construed as related to the stellar mass. In the core- 
accretion mechanism of planet formation, higher disk masses 
can contain an amount of metals available to form giant plan- 
ets even at lower metallicities. 
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